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20.  ABSTRACT  (CONT’D) 

the  set  of  two  partial  differential  equations  is  recapitulated  together  with 
appropriate  boundary  conditions.  For  vibration  problems,  two  sets  of 
eigenvalue  problems  are  formulated  to  satisfy  the  simultaneous  partial 
differential  equations  and  the  homogeneous  boundary  conditions.  Suitable 
parameters  are  defined  to  describe  the  dispersion  relations.  These  dual 
eigenvalue  matrix  equations  are  then  solved  numerically.  For  an  Infinite  rod, 
a  dispersion  relation  of  frequency  versus  wave  number  which  contains  an 
imaginary  branch  has  been  obtained.  The  free  vibration  problem  of  a  fixed- 
fixed  Mindlln-Herrmann  rod  has  been  solved.  The  numerical  values  of  six  (6) 
lowest  frequencies,  the  associated  wave  numbers  and  mode  shapes  are  tabulated 
for  three  different  slenderness  ratios. s 
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INTRODUCTION 


This  report  investigates  the  free  vibration  problem  of  a  rod  of  finite 
length  which  models  the  axial  as  well  as  radial  modes  of  motion.  The  formula¬ 
tion  is  based  on  a  rod  model  of  infinite  length  developed  by  Mindlin  and 
Herrmann  in  1951  (ref  1).  The  differential  equations,  boundary  conditions, 
and  the  associated  energy  principle  have  all  been  developed  in  the  original 
paper  as  well  as  dispersion  relations  (relations  between  wave  number  and 
velocity)  for  a  rod  of  infinite  length.  Subsequently,  Herrmann  presented  a 
solution  formulation  for  free  and  forced  vibration  of  rods  of  finite  length 
(ref  2).  No  numerical  results  were  given  there,  however.  Miklowitz  has 
obtained  formed  solutions  to  forced  vibration  problems  of  rods  of  semi¬ 
infinite  and  finite  length  by  the  use  of  Laplace  transforms  (ref  3).  He  has 
also  obtained  numerical  results  for  the  earlier  formulations  (ref  4).  Never¬ 
theless,  no  specific  numerical  data  were  provided  for  free  vibrations  of  a 
finite  rod.  Since  the  free  vibration  information  is  fundamental  to  all  linear 
and  nonlinear  wave  propagation  phenomena,  our  first  step  is  to  study  this 
aspect  for  a  Mi ndl in-Herrmann  rod  of  finite  length. 


1-R.  D.  Mindlin  and  G.  Herrmann,  "A  One-Dimensional  Theory  of  Compressional 
Waves  in  an  Elastic  Rod,"  Proceedings  of  the  First  U.S.  National  Congress  of 
Applied  Mechanics,  1950,  pp.  187-191. 

^G.  Herrmann,  "Forced  Motions  of  Elastic  Rods,”  Journal  of  Applied  Mechanics, 
September  1954,  pp.  221-224. 

3j.  Miklowitz,  "Travelling  Compressional  Waves  in  an  Elastic  Rod  According  to 
the  More  Exact  One-Dimensional  Theory,"  Proceedings  of  the  Second  U.S. 
National  Congress  of  Applied  Mechanics,  June  1954,  ASME,  1955,  pp.  176-186. 

* J .  Miklowitz,  "The  Propagation  of  Compressional  Waves  in  a  Dispersive  Elastic 
Rod,”  June  1957,  pp.  231-239. 


First,  the  governing  partial  differential  equations  and  pertinent 
variables  of  the  approximate  theory  of  the  rod  are  recapitulated.  With  the 
original  equations,  an  eigenvalue  problem  based  on  a  wave  motion  Is 
formulated.  The  characteristic  equation  of  this  eigenvalue  problem  is  shown 
to  be  identical  to  the  one  in  Reference  l  in  terms  of  phase  velocity  and  wave 
number.  However,  the  same  equation  can  be  written  In  terms  of  dimensionless 
frequency  and  wave  number  which  is  also  allowed  to  assume  imaginary  values. 
Then,  the  dispersion  relation  of  frequency  versus  wave  number  reveals  another 
Imaginary  branch  similar  to  the  one  corresponding  to  Timoshenko  beam  theory 
(ref  5).  This  new  dispersion  relation  is  shown  here  together  with  other 
results  identical  to  the  original  paper  by  Mindlin  and  Herrmann.  Next, 
solutions  for  finite  rods  are  treated.  As  a  first  example,  we  take  a  rod  with 
fixed-fixed  ends.  In  order  to  satisfy  these  end  conditions,  another 
eigenvalue  problem  is  established.  This  produces  a  characteristic  equation 
for  those  particular  frequencies  so  that  free  vibrations  under  the  given 
boundary  conditions  are  possible.  The  numerical  results  are  presented;  and 
among  infinite  number  of  discrete  frequencies,  several  of  the  first  lowest 
are  recorded.  Corresponding  to  each,  the  various  modes  (consisting  of  sine, 
cosine,  hyperbolic  sine,  and  cosine  functions)  and  relative  amplitudes  are 
also  presented  in  the  last  section. 

Ir.  D.  Mindlin  and  G.  Herrmann,  "A  One-Dimensional  Theory  of  Compressional 
Waves  in  an  Elastic  Rod,”  Proceedings  of  the  First  U.S.  National  Congress  of 
Applied  Mechanics,  1950,  pp.  187-191. 

"’S .  P.  Timoshenko,  "On  the  Correction  for  Shear  of  the  Differential  Equation 
for  Transverse  Vibrations  of  Prismatic  Bars,”  Philosophical  Magazine, 

Series  6,  Vol.  41,  1951,  pp.  744-746.  '  "  ~  ' 


GOVERNING  EQUATIONS 


It  will  be  convenient  here  to  recapitulate  the  full  set  of  equations  of 
the  Hindi In-Herrmann  rod  model  (ref  1).  The  displacement  equations  of  motions 
are: 

a2(A+2p)w"  +  2aAu'  +  2aZ  =  pa2w 
a2k2u”  -  8k^2(A+p)u  -  4aki2Aw'  +  4aR  =  pa2u 
where  w  =  w(z,t)  is  the  displacement  component  in  z-direction;  u  =  u(z,t)  so 
that  ru(z,t)/a  is  the  displacement  component  in  r-direction;  z,  r  denote  the 
axial  and  radial  coordinates,  respectively;  t,  the  time.  A  prime  (')  denotes 
partial  differential  with  respect  to  z  and  a  dot  (•)>  the  same  with  respect  to 
t.  A,  p  are  Lame'  elastic  constants,  p  is  the  density  of  the  rod  material, 
and  a  is  the  radius.  R  and  Z  denote  the  radial  and  axial  components, 
respectively,  of  the  traction  on  the  cylindrical  surface  of  the  rod.  k  and  kj 
are  correction  factors.  They  are  introduced  so  that  the  dispersion  relations 
from  this  approximate  theory  will  better  match  the  exact  solution  of  an 
infinite  rod.  Their  usage  is  discussed  fully  in  the  original  paper  by  Mindlln 
and  Herrmann  (ref  1). 

The  stress-displacement  relations  are 

2Pr  =  2P@  *  k^2[2a(A+p)u  +  a2Aw'] 

2PZ  =*  2aAu  +  a2(A+2p)w' 

4Q  =  k2a2pu' 


^R.  D.  Mindlin  and  G.  Herrmann,  "A  One-Dimensional  Theory  of  Compressional 
Waves  in  an  Elastic  Rod,"  Proceedings  of  the  First  U.S.  National  Congress  of 
Applied  Mechanics,  1950,  pp.  187-191. 


where  Pr,  Pq,  Pz,  and  Q  are  the  averaged  stresses  which  are  related  to  the 


Cauchy  stress  components  orr,  etc.,  by  the  following  relations: 


pr  “  J  °rrrdr; 


P0  =  J0  °O0tdr 


(3) 


f3  o 
0 


,a  °rz 


zz 


rdr ;  Q  =  f  ---  r2dr 
0  a 


Equation  (2)  will  be  used  to  obtain  solutions  of  finite  rods  if  stress 
boundary  conditions  at  the  ends  are  given. 

In  the  following  discussion,  we  shall  assume  that  the  cylindrical  surface 
of  the  rod  is  free  of  stress.  Hence 


and  Eq.  (1)  becomes: 


R  =  Z  =  0 

a2(X+2vi)w"  +  2aXu'  =  pa2w 
a2k2u”  -  8ki2(X+u)u  -  4aki2Xw’  *  pa2u 


(4) 


d') 


SOLUTION  FORMULATIONS  FOR  RODS  OF  INFINITE  LENGTH 

To  solve  Eq.  (1'),  we  divide  through  by  pa2  and  obtain: 

X+2u  2X 

( - )w"  +  —  u'  =■  w 

P  pa 

k2u  Ski2(\+\1)  4kx  2X 

( - )u" - y - U - W’  -  U  , 

P  pa  pa 


Let 


iTze-iwt 

w(z,t)  =  Ae 


) 


(3) 


(6) 


i  Yz_-iu>t 


and  substitute  Eq.  (6)  Into  Eq .  (5).  One  has 


X+2  y  ,  2  X 

[-( - )Y  +  ^  ]A  +  1y( — )B  =  0 

P  pa 


4k12X 


k2M  ->  8k]/(X+p) 


-lY( - )  A  +  I—  Yz - 3 - +  u/]B  =  0 

pa  P  pa^ 


Divide  Eq .  (7)  through  by  Y2 : 


( X+2  u)  w2  2X 

[ - 4-  __]A  +  1( - ) B  =  0 


4kt2X  k2p  8k12(X+p)  ^2 

-1( - )  A  +  [ - 5-5—  +  -o]B  =  0 

PYa  P  PY  az  Y 


From  Eq.  (7'),  one  obtains  the  following  characteristic  equation: 

X+2u  w2  k2  p  8k^2(X+p)  u2  8k^2^2 

( - 2)( - f - 2~2 - 2>  +  TH  =  0 

P  Y  P  PY4a  Y  PYa 


and  the  amplitude  ratio 


B  X+2  u  u>2  2  X 

«  =  -  =  -[ - *  -ol /l K — )) 

A  r\  v*-  fivo 


4ki 2  X  k2p  8k^  2(  X+u)  u2 

+  [1( - )]/[ - o'*?  +  ’ll 

PYa  p  PY  a4  Y 


where  Eq .  (8)  Is  exactly  the  same  as  Eq .  (22)  in  Mindlin-Herrmann' s  paper  of 
1951  (ref  1). 

It  would  be  convenient  for  the  present  analysis  to  Introduce  some 
dimensionless  parameters.  But  first,  let 


lR.  D.  Mindlln  and  G.  Herrmann,  "A  One-Diraenslonal  Theory  of  Compresslonal 
Waves  In  an  Elastic  Rod,”  Proceedings  of  the  First  U.S.  National  Congress  of 
Applied  Mechanics,  1950,  pp.  187-191. 


£  =  length  of  the  rod 
L  =  wave  length 


and  hence  the  wave  number  can  be  defined  as 

y  =  2tt/L 

Now,  define 

a*  =  a/£  ,  y*  =  yi  =  2ir£/L 

(w*)2  =  u2a2/(p/p)  ,  =  A/p 

If,  from  Eq.  (6),  we  define  wave  (phase)  velocity  as 

c  =  u/y 

Cg  =  (p/p)1/2 


) 


and 


one  has  the  relation 


w2a2 


P  2  M  2  2 
-  y  -  Y  az 


uj*'1 


Y*  a 


2„*  2 


P  P 

With  Eqs.  (11)  and  (13),  Eq .  (8)  becomes  the  following  when  divided 


P2/p2 


U)* 


8k!2(A*+l)  w*2 


8kL2A*2 


(A*  +  2  -  72“*2^k2  +  ~~~7o  “™)  """  0 


Y*  a 


Y*2a*2  Y*2a*2  Y*2a*2 


Equation  (9)  can  also  be  written  as 

B 


a  -  -  «  -1(A*  +  2 - - >  - 


w*2 

Y*2a* 


2  A* 
Y*a* 


4kL2A*  8ki2(A*+l)  ^2 

/(k2  + - - - - -) 


Y*a 


y*  a* 


Y*^a 


(10) 

(U) 

(12) 

(13) 

through  by 

(14) 

(15) 


Now,  we  shall  solve  Y*a*  for  any  given  uj*  In  Eq.  (14).  Let 


In  terms  of  x  and  y,  Eq.  (14)  becomes 


(A*+2-xy)[k2  +  8k12(A*+I)y  -  xy]  -  8k!2A*2y  =  0 

Or,  upon  expanding, 

[x2  -  8k12(X*+l)x]y2  -  [  ( k2+A*+2 )x  -  8k]_  2(3  A*+2 )  ]  y  +■  k2(A*+2)  =  0 

Or, 

Ay2  -  By  +  C  =  0 

with 

A  =  x2  -  8k12(A*+I)x  ^ 

B  =  (k2+A*+2)x  -  8k12(3A*+2)  , 

C  =  k2( A*+2) 

Thus,  one  can  write  _ 

B  ±  /B2"  -  4AC 


(17) 


(18) 


(19) 


(20) 


Before  we  proceed  to  the  actual  computations,  some  observations  on  the  sign  of 
A,  B,  and  B2  -  4AC  should  be  helpful. 

1.  A  3  x[x  -  8kf 2 ( A*+l ) ] ,  and,  since  x  =  w*2  is  always  positive  or  zero, 
for  w*  *  0,  A  will  be  negative,  zero,  or  positive  depending  on  whether  x  Is 


less,  equal,  or  greater  than  8kf2(A*+l) 


<  <  , 

A  =  0  If  x  =  Ski 2( A*+l) 

>  > 


16 


2. 


Hence 


B  =  (k2+A*+2)x  -  8kf  2(3  A*+2) 


-  (k2+A*+2)[x  - 


8k12(3A*+2) 

k2+A*+2 


<  <  8kL2(3A*+2) 

B  =  0  If  x  - - - - 10 

>  >  k2+A*+2 


where  the  symbol  ~  indicates  "in  the  neighborhood  when  k,  and  X*  are 
recognized  to  have  an  order  of  magnitude  of  unity." 

3.  Finally,  let  us  consider  B2  -  4AC: 

B2  =  (k2+X*+2)2x2  +  64k!4(3X*+2)2  -  16k12(3X*+2)(k2+X*+2)x 
4AC  -  4[x2  -  8k12(X*+l)x]k2(X*+2)  =  4k2(X*+2)x2  -  32k2kt  2(  X*+l )  (  X*+2  )x 
=  4k2(X*+2)x2  -  32k^  2( X*2+3 X*+2)k2 
B2  -  4AC  =  [k2  -  ( X*+2) ] 2x2  +  64kx4(3X*+2)2 
-  16k!2x[3X*+2)(X*+2)  -  k2(2 X*2+3X*+2) 1 

a.  For  x  being  very  large  compared  with  unity,  the  first  term 
dominates,  hence,  B2  -  4AC  >  0. 

b.  For  x  being  very  small  compared  with  unity,  the  second  terra 
dominates,  and  hence,  B2  -  4AC  >  0  again. 

c.  For  x  being  in  the  neighborhood  of  unity,  B2  -  4AC  ~  4  +  1600  -  128  = 
1476  >  0. 

Hence,  in  most  likelihood,  B  -  4AC  >  0  for  all  values  of  x.  But  we  shall 
examine  this  carefully  in  actual  computations. 

NUMERICAL  RESULTS  FOR  RODS  OF  INFINITE  LENGTH 

The  solution  formulations  of  the  previous  section  are  now  carried  out  for 
specific  numerical  parameters.  The  results  are  shown  in  Figures  1,  2,  and  3. 
In  Figure  1,  two  dispersion  curves  (i.e.,  c/c  vs.  a/L)  reported  in  Reference  1 
are  reproduced  here  for  the  purpose  of  numerical  verifications.  As  in 
Reference  l,  we  have  set  v  =  0.29.  The  curve  la  corresponds  to  k^  ■  l  and 

1r.  D.  Mindlin  and  G.  Herrmann,  "A  One-Dimensional  Theory  of  Compressional 
Waves  in  an  Elastic  Rod,"  Proceedings  of  the  First  U.S.  National  Congress  of 
Applied  Mechanics,  1950,  pp.  187-191. 


k  *  0.9258  which  is  obtained  from  Eq.  (25)  in  Reference  1;  the  curve  l'a 
corresponds  to  the  same  k,  but  k^  is  obtained  by  Eq.  (28)  of  the  same 
reference.  In  Figure  2,  the  coupling  effect  of  the  axial  and  radial  modes 
through  Poisson's  ratio  given  in  Reference  1  is  also  reproduced  here.  As 
shown  in  Figure  2,  there  exist  two  velocities  for  each  wave  number  a/L  in  its 
full  range.  However,  it  will  be  realized  that  the  full  range  of  frequency 
(w*)  is  not  covered  in  Figure  2.  Since  <jJ*  is  related  to  the  velocity  c/cs  by 
Eq.  (13),  it  is  necessary  to  solve  u*  for  given  a/L  =  2TfY*a*  in  the 
characteristic  equation  (14).  Hence,  another  dispersion  relation  is  obtained 
and  is  shown  in  Figure  3  for  w*  vs.  2iry*a*  curves.  Evidently,  in  this  figure, 
full  range  of  frequency  is  covered.  It  is  observed  that  for  w*  smaller  than  a 
value  near  4.3,  one  root  of  y*a*  is  imaginary,  while  the  other  is  real.  The 
eigenfunctions  corresponding  to  the  imaginary  root  are  hyperbolic  sine  and 
cosine  functions.  They,  together  with  the  sine  and  cosine  functions 
associated  with  the  real  root  of  Y*a*,  will  be  needed  for  vibration  solutions 
for  a  given  set  of  end  conditions  of  a  finite  rod. 

SOLUTION  FORMULATIONS  OF  RODS  OF  FINITE  LENGTH 

From  solutions  of  rods  of  infinite  length,  one  observes  that  there  are 
four  fundamental  solutions  in  the  form  of  sines  and  cosines  or  hyperbolic  sine 
and  cosine  functions  corresponding  to  any  given  frequency  u>.  Now,  we  wish  to 
find  a  particular  oi  for  a  given  specific  length  of  the  rod  and  for  a  given  set 
of  boundary  conditions.  It  is  also  observed  that  the  fundamental  solutions  of 

1r.  D.  Mindlin  and  G.  Herrmann,  "A  One-Dimensional  Theory  of  Compressional 
Waves  in  an  Elastic  Rod,”  Proceedings  of  the  First  U.S.  National  Congress  of 
Applied  Mechanics,  1950,  pp.  187-191. 
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u)  and  u  are  matched  as  follows: 


where 


H 

3 

cos  Y*z*e“^Uit 

with 

CN 

II  3 

sin  Y*z*e“ia)t 

with 

cosh  Y*z*e-iUit 

with 

II  3 

sinh  Y*z*e”^u)t 

with 

u  =  ui  sin  Y*z*e“iu)t 
u  =  U2  cos 

u  =  113  sinh  Y*z*e-iwt 
u  =  114  cosh  Y*z*e-i(Jt 


(21) 


ui  =  aiwl  •  1  =  i.2.3.4 


(22) 


and 


10*' 


2  X* 


<*1  "  ~[-(X*+2)  +  --“Hi/ (---7) 
Y*za*z  Y*a* 


“2  “  “U1 


10*' 


2X* 


«3  =  [  ( x*+2 )  + 

Y*  za*  z  Y*a* 


«4  =  -a3 


\ 

► 

J 


(23) 


From  the  dispersion  relations  one  observes  that  (Y*a*)2  always  has  two  roots 
for  any  given  to*:  either  (1)  both  of  these  roots  are  real,  or,  (2)  one  is 
real  and  the  other,  imaginary. 

For  case  (1),  the  general  solution  has  the  form 

3  3  =  * 

w(z,t)  *  (wi  cos  Yi*z*  +  W2  sin  Yi*z*  +  W3  cosh  Y*z*  +  W4  sinh  Y*z*)e"*U)t  (24) 

3  3  3  3 

u(z,t)  =  (ui  sin  Yi*z*  +  U2  cos  Yi*z*  +  U3  sinh  Y*z*  +  U4  cosh  Y*z*)e”*-Wt  (25) 

As  our  first  example,  the  boundary  conditions  of  a  rod  with  both  ends 
fixed  will  be  considered: 


For  a  nontrivial  solution  of  Eq .  (27),  one  has 


Or, 


1 

0 

1 

0 

cos  Y^* 

sln  Yi* 

cos  Y3* 

sln  Y3* 

0 

a2 

0 

a4 

a3  sin  Y]_* 

a2  cos  Yi* 

«3  sin  Y3* 

«4  cos  Y3' 

sin  Yi* 

cos  Y3*  - 

cos  Yi* 

sln  Y3* 

“2 

0 

04 

a2  cos  Yi* 

013  sin  Y3* 

~  <*1  sin  Yi* 

04  cos  Y31 

=  a2  sin  Y3*(«3  sin  Y3*  -  ax  sin  Yi*) 
+  “2a4  sln  Yl*(“3  sin  Y3*  -  sln  Yi*) 

-  «4  sln  Yi*(ct3  sln  Y3*  -  sin  Yi*) 

-  a2a^  cos  Y3*(cos  Y3*  -  cos  Yi*) 


A1  =  (a3  sln  Y3*  -  sin  Yi*)(a2  sln  Y3*  -  %  sln  Yi*) 
-  cvja^Ccos  Y3*  -  cos  Yi*)2  =*  0 
For  a  nontrivial  solution  of  Eq.  (29): 


— 

— 

1 

0 

1 

0 

A  2  * 

cos  Yi* 

0 

sln  Yi* 

a2 

cosh  Y3* 

0 

slnh  Y3* 

04 

a3  sin  Yi* 

a2  cos  Yl* 

»3  sinh  Y2* 

04  cosh  Y3* 

m 


k.:-- 


m 


V  ' 


I  > 


% 


[<V 


F «N  ■ 

£? 


/,  , 

>**J 


•s 
»*  i* 


V*J\ 


Or, 


^2  ■  (013  si"h  Y3*  ~  sin  Yi*)(a2  sinh  Y3*  -  C14  sin  Y3*) 

-  o^a^cosh  Y3*  -  cos  Y3*)2  =  0  (34) 

With  Eqs.  (33)  and  (34),  one  can  write 


A^(oj*)  if  Y*2a*2  has  two  real  roots 


A  =*  A(uj* ) 


A2(w*)  if  Y*2a*2  has  one  real  and  one  imaginary  root 


(35) 


>  O 

Hence,  for  a  given  w*,  one  must  first  find  two  values  for  Y*  a*  and  u^,  a2, 
«3,  and  04  (or  03  and  04).  Then  with  given  a*  =  a /£,  one  obtains  Yi*  and  Y3*. 
Thus  one  can  plot  (w*)  with  respect  to  oj*  and  find  w*’ s  at  which  A(u*)  =*  0. 

Finally  in  this  section,  we  will  write  down  the  formulas  for  computing 
the  relative  amplitudes  of  various  mode  shapes.  Equation  (27)  for  case  (1) 
can  be  written  as 


0 

sin  Y]_* 

1 

cos  Y3* 

0 

sin  Y3* 

w2 

=s 

W3 

1 

COS  Y3* 

a2 

0 

«4 

a 

w4 

0 

w 


(36) 


From  Eq.  (36),  one  has 


D2 


w2  “  w^  =  32w^ 


d3 


W3  »  —  W3  =  33W1 


(37) 


D4 


w4  “  “  wi  -  ^AWl  / 
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*  4  -» 


where 


0 

1 

0 

D  = 

sin 

Yi* 

cos  Y3* 

sin 

"*3* 

=  02  s  in 

Y3*  -  04  sin  Yi* 

a2 

0 

“4 

— 

-1 

1 

0 

— 

D2  =* 

-cos 

1  Yi* 

COS  Y3* 

sin 

y3* 

=  a^(cos 

Yi*  -  cos  Y3*) 

0 

-1 

a4 

— 

0 

-1 

0 

D3  = 

sin 

Yl* 

-cos  Yi* 

sin 

Y3* 

=  04  sin 

Yi*  -  02  sin  Y3* 

u2 

0 

0 

— 

0 

1 

-1 

— 

D4  = 

sin 

Yi* 

COS  Y3* 

-cos 

Yi* 

*  02 (cos 

Y3*  -  cos  Yi*) 

“2 

0 

0 

(38) 


For  case  (2),  the  above  equation  can  still  be  used  with  cos  Y3*  and  sin 

3 

Y3*  replaced  by  cosh  Y3*  and  stnh  Y3*.  Now,  for  u^,  i  =  l, 2, 3, 4,  we  observe 
Eqs.  (22),  (23),  (28),  and  (30),  and  write 

3  3  —3 

ui  3  uiwi  3  &iwi  »  1  3  1.2, 3, 4  (39) 

with 

3i  3  ai3i  ,  i  =  1,2, 3, 4  (40) 
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dalAli 


* .  v  •  •  V  *  »  ' «  *  ^ 


where  3^  =  1  when  w^  is  used  as  the  basis  of  normalization, 

»i 

Pi  =  —  ,  i  =  2,3,4 

Dl 

(41) 

h  =  1 

In  Eq.  (40),  03  and  04  will  be  used  in  place  of  03  and  04  if  the  mode  shape 
involves  hyperbolic  cosine  and  sine  as  stated  in  the  earlier  discussion. 


NUMERICAL  RESULTS  ON  VIBRATIONS  OF  FINITE  RODS  WITH  AXIAL  AND  RADIAL  MOTIONS 
The  formulations  obtained  in  the  previous  section  are  now  used  for 
specific  numerical  examples.  In  conjunction  with  a  fixed-fixed  rod  (i.e.,  w 
and  u  are  zero  at  both  ends),  several  lowest  vibration  frequencies  and  the 
associated  eigenfunctions  have  been  obtained.  The  results  are  tabulated  in 
Tables  I  through  III  for  slenderness  ratio  a*  =  a/ &  =  0.05,  0.10,  and  0.20, 
respectively.  In  these  tables,  3^,  i  =*  1,2, 3, 4  are  the  coefficients  for 
w(x,t)  and  Bj  for  u(x,t).  They  are  normalized  with  respect  to  the  one  with 
largest  absolute  value.  For  any  given  u)*,  a  pair  of  Y*a*  is  obtained.  These 
are  indicated  by  Y\*a  and  Y3*a,  respectively,  as  described  in  the  previous 
section.  If  one  of  them  (Yi*a)  is  imaginary,  i.e.,  Yi*a  =  iYi*a,  then  Yi*a  is 
given  and  it  is  indicated  by  parentheses.  The  coefficients  of  the  associated 
eigenfunctions  (hyperbolic  sine  and  cosine  functions)  are  also  Indicated  by 
parentheses.  For  example,  in  Table  I  for  a*  =  a/ £  =  0.5,  both  the  lowest 
frequencies  have  one  imaginary  root  each.  For  all  higher  aJ*’s,  both  roots  of 
Y*a*  are  real. 
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VIBRATION  SOLUTIONS  TO  A  FINITE  MI NDL IN-HERRMANN  ROD  WITH  FIXED-FIXED  END  CONDITIONS 

(v  *  0.29,  Slenderness  Ratio  =  0.05) 
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The  numerical  values  in  parentheses  are  associated  with  imaginary  wave  numbers  and  their 
hyperbolic  sine  and  hyperbolic  cosine  (mode  shape)  functions. 


VIBRATION  SOLUTIONS  TO  A  FINITE  MINDLIN-HERRMANN  ROD  WITH  FIXED-FIXED  END  CONDITIONS 

(v  -  0.29,  Slenderness  Ratio  **  0.10) 


NOTE:  The  numerical  values  in  parentheses  are  associated  with  imaginary  wave  numbers  and  their 
hyperbolic  sine  and  hyperbolic  cosine  (mode  shape)  functions. 


VIBRATION  SOLUTIONS  TO  A  FINITE  MI NDL IN-HERRMANN  ROD  WITH  FIXED-FIXED  END  CONDITIONS 
(v  -  0.29,  Slenderness  Ratio  =  0.20) 
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The  numerical  values  in  parentheses  are  associated  with  imaginary  wave  numbers  and  their 
hyperbolic  sine  and  hyperbolic  cosine  (mode  shape)  functions. 


Mindl in-Herrmann  theory  with  <  correction 
Mindl in-Herrmann  theory  with  <  and  k. 


Longitudinal  wave  velocities  of  first  mode  of  motion  according 
to  Mindl in -Herrmann  Rod  Theory  (same  as  curves  la  and  l'a  in 
Figure  1,  Reference  1). 


2.  Effect  of  coupling  through  Poisson's  Rati 
wave  velocities  of  axial  and  radial  modes 
(same  as  Figure  2,  Reference  1). 


Im(2?rr  a  ) 


Figure  3 


Frequency  versus  real  and  imaginary  "wave  numbe 
in  an  infinite  Mindl in-Herrmann  Rod. 
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